Quantum Orders and Spin Liquids in Cs 2 CuCl4 
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Motivated by experiments on CS2CUCI4 samples, we studied and classified the symmetric spin 
liquids on triangular lattice. We identified 63 Z2 spin liquids, 30 U(l) spin liquids and 2 5(7(2) spin 
liquids. All those spin liquids have the same symmetry but different quantum orders. We calculated 
the spin spectral functions in some simple spin liquids and compared them to the one measured 
by experiments on CS2CUCI4. We find that the UICt+t'Lt 1 spin liquid or one of its relatives is 
consistent with observed properties of the spin liquid state in CS2CUCI4. We discussed the fine 
distinctions among those possible spin liquids and the spin liquids proposed using slave-fermion 
approach, so that future experiments can determine which spin liquids actually describe the spin 
state in CS2CUCI4. 

PACS numbers: 73.43.Nq, 74.25.-q, 11.15.Ex 
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I. INTRODUCTION 

Our understanding of states of matter and their inter- 
nal orders has been dominated by Landau's symmetry 
breaking theory. 0, Q We used to believe the all possi- 
ble orders are described by various symmetry breaking 
states. However, after the discovery of fractional quan- 
tum Hall (FQH) states, H || we started to realize that 
FQH states contain a new kind of order - topological 
order 5] (which was first proposed to describe spin liquids 
found in research of high T c superconductors^). What 
is new about topological orders is that topological or- 
ders cannot be characterized by symmetry breaking and 
the related local order parameters and long range cor- 
relations. Topological orders are characterized by new 
set of universal quantum numbers, such as ground state 
degeneracy, fractional statistics, edge excitations, eic.[jj 

Recently, a concept of quantum order was proposed^ 
8] to describe non-symmetry breaking order that gener- 
ally appear in a quantum state. The quantum order gen- 
eralize the topological order to gapless states. To have a 
concrete description of quantum order without using or- 
der parameters, a new mathematical object - projective 
symmetry group (PSG) - was introduced. The concept 
of quantum order and its PSG characterization allow us 
to understand quantum phases and quantum phase tran- 
sitions in a systematic way. 

Just like the group theory allows us to classify sym- 
metry breaking orders, the PSG characterization of 
quantum states allows us to classify different quan- 
tum orders^ and can distinguish two different quantum 
phases even when they have exactly the same symmetry. 

Also just like the symmetry description of classical or- 
der allows us to obtain low energy properties of system 
without knowing the details of the systems, the PSG de- 
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scription of quantum orders also allows us to obtain low 
energy properties of system without knowing the details 
of the systems. However, unlike symmetries whichpro- 
duce and protect gapless Nambu-Goldstone bosons, 00] 
the quantum orders can produce and protect gapless col- 
lective modes which behave like light and other gapless 
gauge bosons. [II], [13] Those gapless collective modes can 
also be gapless/massless fermions, even when the origi- 
nal theory is purely bosonic.[l3l[l4],[l5|] More recently, it 
was realized that the quantum ordered states described 
by the PSG are actually string- net condensed states. [T^ 
The emerging gauge bosons are the fluctuations of con- 
densed string-nets and emerging fermions are the end of 
condensed strings. [T^ 

In this paper, we would like to further develop the PSG 
description of quantum orders. In Rcf. 7], using PSG, a 
classification of quantum orders in symmetric spin liq- 
uids on 2D square lattice was given. Here, we would like 
expand the results of Ref.Q to 2D spin liquids which 
break the parities P x : x — > —x and P y : y — > —y. We 
will assume the spin liquids have the following symme- 
tries: two translation symmetries T x , T y , two parities 
P xy : (x,y) -> (y,x) and P xy : (x,y) -> (-y,-x), and a 
time reversal symmetry. Such type of 2 D sp in liquids was 
observed recently in CS2CUCI4 sample, [la, Il9j By classi- 
fying the quantum orders in those spin liquids, we hope 
to identify the quantum order in the CS2CUCI4 sample. 

The spin dynamics in CS2CUCI4 can be described by 
a 2D spin- 1/2 system on a square lattice with nearest 
neighbor coupling J' = 0.125meV and one diagonal cou- 
pling J = 0.375meV in the x + y direction. At temper- 
ature T > T c — 0.62K= 0.053meV, the 2D spin system 
was found to be in a liquid state. Since T c is a small 
energy scale, we will regard the finite temperature spin 
liquid state as a zero temperature quantum state. (More 
precisely, we assume that we can add additional frus- 
trations to lower T c to zero.) The theoretical spin-wave 
calculations [23, [2l| and series expansion calculations [22j 
also suggest the existence of spin liquids in the above J- 
J' model. In this paper, we are going concentrate on the 
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following physical issue: what is this spin liquid in the 
CS2CUCI4 sample. There are several possibilities. 

Since J' is small, one possibility is that the spin liq- 
uid, at low energies, just behaves like a decoupled ID 
spin-1/2 chains. In other words, the spin liquids is in 
the same universality class of the decoupled ID spin-1/2 
chains. It is also possible the spin liquid is not purely 
ID and has intrinsic 2D correlations. The 2D correla- 
tion can lead to 2D dispersion for low lying excitations. 
Furthermore, there can be several different types of 2D 
spin liquids that have the same symmetry as the purely 
ID spin liquid. Since all those possible spin liquids have 
the same symmetry, it is difficult to study them without 
knowing how to cgaracterize them. 

In this paper, using quantum order and its PSG, we 
construct and characterize a large class of 2D spin liquids 
that have the same symmetry as the purely ID spin liq- 
uid. In section^ we introduce mean-field ansatz that de- 
scribe the symmetric spin liquids - the spin liquids that do 
not break any symmetries. In section ITTT1 we discuss how 
to use PSG's to characterize different mean-field phases 
(or the universal classes of the mean-field ansatz). In 
section |IV| we find all the Z 2 and SU(2) PSG's and a 
large class of U(l) PSG's within the SU(2) slave-boson 
approach. This tells us the possible spin liquids that can 
be constructed use the SU(2) slave-boson theory. 

Constructing a large class of symmetric spin liquids 
and obtaining their PSG characterization are useful in 
the following sense. If a spin liquid state is found in 
CS2CUCI4 or some other samples which do not break 
any symmetry, then the spin liquid has a good chance 
to be in the class that we obtained. Identifying the PSG 
that characterizes the constructed spin liquid will allow 
us to identify many universal properties of the spin liquid. 
We can check those universal properties experimentally 
which allow us to identify the PSGs for the experimen- 
tally observed spin liquids. The classification will also 
help us the study the phase transitions between symmet- 
ric spin liquids that do not change any symmetries. 

In section ]S\ we discuss some simple ansatz that re- 
alize some of the classified mean-field spin liquids. We 
calculate the mean-field energies of the constructed spin 
liquids to determine which spin liquids are likely to ap- 
pear in the CS2CUCI4 system. We also calculate spin cor- 
relation in those likely spin liquids in section lVll This al- 
lows experimentists to determine which spin liquid actu- 
ally describe CS2CUCI4 system (above T c ) using neutron 
scattering and other techniques. Section IVIII discusses 
connection of our results with experiments and previous 
theoretical results. 

The main point of our calculation is to identify univer- 
sal properties of various spin liquid states. This point is 
highly non-trivial since the involved spin liquids all have 
the same symmetry. We show that the crystal momenta 
of gapless spin-1 excitations can be used to experimen- 
tally distinguish different spin liquids. 

We would like to pointed out that there are two ways to 
study possible spin liquid states in CS2CUCI4. The first 



approach is the slave-fermion approach, where one start 
with a spin ordered state (such as the spiral state) and 
then study the spin liquid induced by strong spin wave 
fluctuations. [23LI24L I25I One can obtain the spin spectral 
function from the slave-fermion approach. The calcu- 
lated spin spectral function 124. [251 agrees well with ob- 
served spin spectral function. |l8l. 1 19| In the slave-fermion 
approach, the spin disordered state (the spin liquid state) 
always has a gap, with low lying bosonic excitations. 

In this paper, we are going to use slave-boson approach 
to study spin liquid states. In the slave-boson approach, 
the spin liquid state can either be gapped or gapless. 
The low lying excitations are fermions. In general slave- 
boson approach can generate more exotic states then the 
slave-fermion approach. We will see later that the spin 
spectral function obtained from the slave-boson approach 
also agrees well with observed spin spectral function. [IH 

m 

However, the spin liquid obtained in Ref. jU 0] and 
the spin liquids obtained in this paper are really different. 
For one thing, the spin liquid obtained in Ref. [H [24| 
has a gap, while the spin liquids obtained in this paper 
is gapless. Because the gap is small, the two spin liquids 
behave similarly at finite temperatures. 

What is the relation between the slave-boson approach 
and the slave-fermion approach? To understand the re- 
lation, we need use the result in Ref. fl6| where it was 
shown that the spin liquids obtained from the slave-boson 
approach and the slave-fermion approach have a conden- 
sation of nets of closed strings. The quantum orders in 
the spin liquids obtained by both slave-fermion and slave- 
boson approaches can be described by PSG's. PSG's 
simply characterize different string-net condensations. 

In the spin liquids obtained by slave-fermion approach, 
we concentrate on one type of condensed strings whose 
ends are bosons (we will call those strings bosonic 
strings). The PSG is nothing but the symmetry group 
of the hopping Hamiltonian of the ends of the condensed 
strings, [lq For different string condensations, the PSG 
for the ends of the condensed strings are different. Hence 
PSG can be used to characterize different string con- 
densations (or different quantum orders), [lq 1 In slave- 
fermion approach, we concentrate on spin liquids with 
condensation of fermionic strings. The ends of fcrmionic 
strings are fermions. We can use the PSG for the ends 
of fermionic strings to characterize different string con- 
densed (or quantum ordered) states. 

Some spin liquids contain condensation of both bosonic 
strings and fermionic strings. Those states can be 
characterized by either the PSG for the ends of the 
bosonic strings or the PSG for the ends of the fermionic 
strings. 16] Those states can be constructed either by 
slave-fermion approach or slave-boson approach. 
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FIG. 1: The nearest neighbor spin model on a triangular lat- 
tice. 



II. SYMMETRIC SPIN LIQUIDS ON A 
SPIN-1/2 SYSTEM 

We consider the spin- 1/2 system on a 2D square lat- 
tice with nearest-neighbor coupling J' and next-nearest- 
neighbor coupling J on only one of the diagonal link 
(ix,iy) — * [ix + + !)■ Such a lattice can also be 
regarded as a triangular lattice which does not has the 
60° rotation symmetry. Within the SU(2) slave-boson 
approach [2^ (23, general spin wave function can be con- 
structed by introducing a mean-field Hamiltonian 

Hmean = ~ ^ (^juijl/jj + h.c)j + ^ a oV4 ^Ipi (1) 
(i,3) 1 

where ip T = (ipi 7 ip2), u]i — Uij, r 1,2,3 are the Pauli ma- 
trices, and Uij are 2x2 complex matrices. The collection 
(uij,a l T 1 ^ is called a mean-field ansatz For each 

mean-field ansatz, we can obtain a mean-field ground 
state by filling the negative energy levels of H mean with 

the spinons tp : \^mean° T } ■ The physical spin wave 
function ^f S pin can now be obtained by performing a pro- 
jection (see 0): 

^^({ri}) = (0| Ui^^r^ti&'h (2) 

i=l 

where J-j is the coordinate of the i th up-spin and N up 
is the total number of the up-spin. Here an empty site 
of ip represent a down-spin and a double occupied site 
represent a up-spin. Single occupied sties are unphysical 
states and are projected out. 

Our representation does not have explicitly spin rota- 
tion symmetry and it is hard to see spin rotation symme- 
try from the mean-field ansatz. This problem is solved 
in Ref. Q. We find that a spin liquid is spin-rotation 
symmetric iff tiy satisfies 

= ipijWij 

Pij = non-negative real number (3) 
W i:i G 577(2) 



In this paper, we will only consider spin liquids with spin 
rotation symmetry, and we will assume that the ansatz 
satisfy the above conditions. 

If we regard the ansatz (Uij , a l T l ) as a label of physical 

(■ui S , a" t 1 ) 

wave function, ^p^' ° then such a label is not 

a one-to-one label. It has a SU(2) gauge structure, i.e. 
two ansatz, (uij,a T l ) and (W(uij),W(a l T 1 )), related 
by a SU(2) gauge transformation label the same physical 
wave function: 

* spm ({aJ) = (0|n/*a 4 |*& ),W(o, ° T,)) > 

i 

= (oin/^i*— ! ° r!) > w 

i 

where W{u l3 ) = W %Uij w] , W(a l Q (i)T l ) = W t a l {i)T l w\ , 
and Wi G SU(2). 



III. QUANTUM ORDER AND PROJECTIVE 
SYMMETRY GROUP 

For the symmetry breaking states, the symmetry is a 
universal property shared by all the states in the same 
phase. Because of this, the symmetry provides a quan- 
tum number that characterize different symmetry break- 
ing orders. Hence we can say symmetry group (SG) de- 
scribes the internal orders of symmetry breaking states. 
To characterize and classify different quantum orders 
which contain no symmetry breaking, we should con- 
struct some universal quantum numbers which describe 
different classes of quantum entanglement in the many- 
body ground state wave function. Ref. and @ propose 
that the symmetry of the mean-field ansatz (uij,a l Q T 1 ^ is 
a universal property and serves as a quantum number 
that characterize the quantum order in spin liquids. The 
symmetry of group of the ansatz will be called the pro- 
jective symmetry group (PSG). An element of PSG is a 
combined operation of a symmetry transformation fol- 
lowed by a gauge transformation. By definition, a PSG 
is formed by all the combined operations that leaves the 
ansatz unchanged. 

Because of the 577(2) gauge structure, the (v,ij,a l T l ) 
labeling of the physical spin wave function 4" S pin({&i}) 
is not a one-to-one labeling. Two mean-field ansatz dif- 
fered by a SU(2) gauge transformation give rise to the 
same spin wave function. Because of this, it is a non- 
trivial task to find out the symmetry of a spin liquid 
label by an ansatz (uij,a l T l ). In order for a spin wave 
function to have a symmetry, its corresponding ansatz is 
only required to be invariant under the symmetric trans- 
formation followed by a proper SU (2) gauge transforma- 
tion. Thus given two spin wave functions with the same 
symmetry, their ansatz can be invariant under the same 
symmetric transformations followed by different gauge 
transformations. In this case, the two spin liquids with 
the same symmetry can have different PSG's. We see 
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that the PSG characterization is more refined than the 
symmetry group characterization. PSG contains some 
information about the phase of spin wave function. 

To understand the precise relation between PSG and 
SG, we need to introduce a subgroup of PSG - the invari- 
ant gauge group (IGG).0j3] IGG is a special subgroup 
of PSG, which is formed by pure gauge transformations 
that leave the ansatz unchanged, 

IGG = {WilWiUijW] = uy, Wt G SU(2)} (5) 

Then SG and PSG are related by SG = PSG/ IGG. 

Spin liquids may support a special kind of low en- 
ergy collective excitations - gauge fluctuations. (2^ It was 
shown thatQ, the gauge group of those low energy 
gauge fluctuations is nothing but the IGG of the corre- 
sponding ansatz. We see that IGG of an ansatz is very 
important. In this paper, we will consider only three 
kinds of IGG, SU(2), U(l) and Z 2 . Correspondingly, we 
call the corresponding spin liquids SU(2), U(l) and Z 2 
spin liquids. 

Mathematically, PSG is defined as, 

PSG = {Gu\GuU{u l3 ) = Uij ,Gu(i) G SU(2)} (6) 

where J7(uy) = fiy = uu(i),uu) > GuU(u tj ) = 
Gij{i)uijG\ J {j) 1 U generates the symmetry transforma- 
tion and Gjj is the associated gauge transformation. 

A PSG may change under a gauge transformation W . 
From WGu(uij) = W( Ulj ), where W{ Uij ) = W iUij w] , 
we find that WG u UW- 1 W(u ij ) = W(u«). Therefore if 
GuU is in the PSG of ansatz uy, then \VGuUW~ 1 is 
in the PSG of the gauge transformed ansatz W(uij). We 
see that the gauge transformation Gjj associated with 
the transformation U changes in the following way 

Gu{i) W(i)Gu(i)W{U(i))l (7) 

under a SU (2) gauge transformation. 

Since PSG is a property of an ansatz, we can group all 
ansatz sharing the same PSG together to form a class. 
Such a class is the universal class of quantum states that 
corresponds a quantum phase. 

Now let us consider some simple examples of spin liq- 
uids described by the ansatz (v,ij,a T J. The first exam- 
ple is 

u i>i+m = u m = u l m T l (8) 

where u l m are real, I = 1, 2, 3. It is easy to obtain the 
spinon dispersion of such an ansatz, 



E ± (k)=±^<-<) 2 

where u l k = J2m u l m e tk ' m . The Brillouin zone is k x ,k y G 
(-%,%). 



The second example is 

Wi,i+x = *XT° + r/r 3 , 

Ui, i+ y = (-)*" (*XT° + vr 3 ) , (9) 

Ui,i+x+y = (—) " ^T 1 , 

where x, rj and A are real. Its spinon spectrum is deter- 
mined by 

H(k) = (x sin k x r° + 77 cos k x T 3 ) t 1 

+ (x sin k y T° + 77 cos k y T 3 ) (g> r 3 
+XT 1 ® r 2 sm(k x + k y ) + a l Q T l (g> r° 

where the Brillouin zone is k x e (—7r/2,ir/2),k y G 
(— 7r,7r). The four bands of spinon dispersion have a 
form of ±Ei(k),±E 2 {k). It seems strange to find that 
the spinon spectrum is defined only on half of the lattice 
Brillouin zone. However, this is not inconsistent with the 
translation symmetry since the single spinon excitation is 
not physical. Only two-spinon excitations correspond to 
physical excitations and their spectrum should be defined 
on the full Brillouin zone. The two-spinon spectrum de- 
fined on the full Brillouin zone can be constructed form 
single-spinon spectrumQ 

E 2s (k) = E ai (h) + E a2 (k 2 ) 

k = ki + k 2 + nirx (10) 

where n = 0, 1 and a± i2 = 1,2 is the sub-index of the 
single spinon dispersion ±E\{k), ±E 2 (k). We note that 
the physical spin-1 excitations are formed by two-spinon 
excitations in the mean-field theory. Thus the The two- 
spinon spectrum E% s is also the spectrum of physical 
spin-1 excitation and can be measured in experiments. 



IV. CLASSIFICATION OF SYMMETRIC SPIN 
LIQUIDS 

In this section, we classify symmetric spin liquids 
on triangular lattices. We consider the spin liquids 
which are invariant under translation transformation 
T x (i — > i + x) and T y (i —> i + y), parity transformation 

Pxy ((^x^ iy) * (^y?^)) &nd P X y{^{i x ^iy) > ( iy? ^x))) 
spin rotation transformation, spin-parity transformation 
T* (S x — ► S x ,S y — > —S y ,S z — ► S z ). The spin rotation 
symmetry requires the ansatz take the form of Eq.©. 
Under the spin-parity transformation T* , the three com- 
ponents of spin change as S x — > S x , S y — > — S y , S z — » 
S^. This transformation can be realized through wave- 
function transformation $ — ► $*, 

under the 5^ representation / and f' are real, then 
(S) = 0y = ~ J <P*paf$ 
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\ UpafQ* = Q / <I>./^/<I> 



which leads to — > S^, 5-, 



ansatz wy transfer as ity 
transformation u.y 
combined transformation as 



Sy, S Z 



S z . The 
combined with a gauge 
(it 2 ) Uij (—ir 2 ), we denote this 



ij ■ 



(11) 



Let us firstly discuss the method to classify the PSG 
briefly and leave the details to the appendix. For two 
given symmetry transformations, whose corresponding 
elements in a PSG, must satisfy some algebraic relations. 
Solving these equations allows us to construct a PSG 
which will be called the algebraic PSG. The new name 
algebraic PSG is introduced to distinguish them from the 
invariant PSG defined in the previous section. Any in- 
variant PSG will be algebraic PSG. But an algebraic PSG 
may not be an invariant PSG unless there does exist an 
ansatz such that the algebraic PSG is the total symmetry 
group of the ansatz. 

As an example, let us consider the two translations T x 
and T yi which satisfy the following relation 



TxTyT-%; 1 



1. 



(12) 



From the definition of PSG, we find that the two elements 
of PSG, G X T X and G y T y: must satisfy 



G X T x Gy Ty {G X T x ) ( Gy Ty ) 

Grp s~l rp rp — \s~i—\rp — \ — 1 
x x^-Xy -*■ y -*- x x y y 

G x (i) Gy (t - x) G- 1 (i - y) G^ 1 (t) e Q. (13) 



Hereafter we will denote the IGG as Q. Each solution of 
the equation corresponds to an algebraic PSG for 
T x and T y . 

By adding other symmetry transformations, we can 
find and classify all the algebraic PSG associated with a 
given symmetry group. Since an invariant PSG is always 
an algebraic PSG, we can check whether the algebraic 
PSG is an invariant PSG through constructing an explicit 
ansatz mj. If an algebraic PSG supports an ansatz «,j 
with no addition symmetry, then it is an invariant PSG. 
Through this method, we can classify the symmetric spin 
liquids through PSG. 

In the appendix, we classify spin liquids with the above 
mentioned symmetries through PSG's. We limit our- 
selves to SU(2), U(l) and Z 2 spin liquids (we only con- 
sider those PSG's whose IGG is one of SU(2), U(l) and 
Z2 )• We found that there are 63 kinds of Z2 spin liquids, 
30 kinds of U(l) spin liquids and 2 kinds of SU(2) spin 



liquids with T XiV , 
tries. 



Pr 



T* and spin rotation symme- 



V. APPLICATION TO NEAREST-NEIGHBOR 
SPIN COUPLING MODEL 

In this section, we present our classification results 
via some simple examples. We will assume only Uij+x, 
ti^i+xj and Ui t i+ x +y are n0 zero. In this case, we have the 
following 7 kinds of Z2 spin liquids and 3 kinds of 17(1) 
(Other Z2 and U(l) spin liquids require non- vanishing 
on longer bonds.) 

The 7 Z2 spin liquids are: 
(1) Z2At {) t\t\ spin liquid 



U lti+X = + VT 2 , 

u i,i+y = XT 1 +VT 2 , 



y — At 1 , 



ai,a = a 2 ,a = 



(2) Z2At x t\t\ spin liquid 



Ui,t +X = + VT 2 , 

Ui,i+v = XT 1 ~ VT 2 , 



= At 1 , 



1 2.3 n 

a n = 01, a = U 



(3) Z2AT Tlrl spin liquid 



Ui,i+£ = ixT° + ?/t 3 , 

Ul,i+y = iXT° + VT 3 , 



(4) Z2At 1 t 3 _t 3 _ spin liquid 



At 1 



2,3 n 

ai,a =0 



,.; + .t = ixT° +VT S , 



u+S = i-XT° - VT 3 , 



a 
u 

^i,i+x+y 



= At 1 , 



1 2.3 n 

a n = 01, a = U 



(5) Z2BT A T n _T\ spin liquid 



Ui,i+ X = X^ + VT 2 , 
Uui+y = (-f* {-XT 1 - VT 2 
Ui,i+x+y = (—) X At 1 , 



1,2,3 








(6) Z2Bt 1 t±t1 spin liquid 



Ui,i+x = XT 1 + VT 2 , 

Ui,i+y = (-f* (xt 1 - VT 2 ) 

Ui,i+x+y = ( — ) X At 2 , 

1.2,3 n 

a = 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 
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(7) ZIBr'Ti.Tt spin liquid 



(20) 



UiA+x+y = (— Y" At 1 , 



1,2,3 



0. 



The 3 U(l) spin liquids are: 
(1) U1At°tItI spin liquids: 





= XT 3 , 




= XT 3 , 


i,i+x-\-y 


= At 3 , 


1,2 

a 


= 0,a 3 = 



(21) 



(2) UIBt 1 t°_t\ spin liquids: 



Ut,l+x = XT 3 , 
Ui,i+y = -(-) l!C XT 3 , 

Ui i+x+t, = (— )* x At 3 , 



(22) 



1,2,3 
V 



= 0. 



(3) U\Ct%t°t 1 spin liquids: 





= x- 2 , 




= -XT", 




= At 3 , 


1.2 

a 


= 0,a 3 = 



(23) 



a 3 . 



The labeling scheme of these spin liquids is defined in the 
appendix. 

We have performed the self-consistent mean-field cal- 
culation for our J-J' model. We found many self- 
consistent mean- field solutions. The energies of four of 
them are plotted in Fig. 2, where we have set J + J' = 1. 
To understand the physical properties of those mean-field 
states, we need to find their PSGs. 

The A phase in Fig. 2 is a one-dimensional state. Its 
IGG is SU oc (2), one SU(2) for each decoupled ID chain. 
We will call such a state SU°°(2) spin liquid. Its ansatz 
is given by 



At ,a 



1,2,3 



u i,i-\-x-\-y 

with spinon dispersion 

E(k) = ±2A cos (k x + k y 



0. 



(24) 



(25) 



This phase, favored by small J', corresponds to the phase 
of decoupled ID spin chain mentioned earlier. 

The B phase is the UIBt 1 ^^ spin liquid. Its four 
spinon bands are given by 

± 2yx 2 (cos 2 k x + cos 2 k y ) + A 2 sin 2 (k x + ky) (26) 




FIG. 2: The mean-field energies for various phases in the 
nearest neighbor spin coupling system. (A) the SU(2) x SU (2) 
spin liquids in Eq.J2]}. (B) the UIB^tItX state in Eq.JSSJ. 
(C) the UICtItIt 1 state in Eq.JUJ. (D) the vr-flux phase 
(SU2B state) in Eq.J^EJ. 




FIG. 3: (a) Contour plot of the spinon dispersion E+(k) as 
a function of (k x /ir,k y /n) for the UIBt^t^t^. state, (b) 
Lower edge of the two-spinon spectrum E2 S (k) as a function 
of (k x /-K,k y /Tx) for the U1Bt 1 t^t\ state. 



The spinon dispersion mm(Ei(k), E2(k)) and the lower 
edge of the spectrum E2 S (k) of the physical spin-1 exci- 
tations are plotted as functions of (kx/ir, k y /ir) in Fig. 3. 
This phase is one of 2D spin liquid phases. 

The C phase has a PSG UICtIt^t 1 . The spinon dis- 
persion is given by 

E + (k) = 2y / [Acos (k x + k y ) — a^] 2 + x 2 (cosk x — cosk y ) 2 

(27) 

The spinon dispersion E + (k) and the lower edge of 
the spin-1 spectrum E2 S (k) are plotted as functions of 
(kx/ir, ky/ir) in Fig. 4. This phase is another 2D spin 
liquid phase. 

The D phase in Fig. 2 is the 7r-flux phase (SU2B spin 



(a) 



(b) 



(a) 



(b) 



FIG. 4: (a) Contour plot of the spinon dispersion E+(k) as 
a function of (k x /n, ky/n) for the UICt+t^t 1 state, (b) 
Lower edge of the two-spinon spectrum E2 S (k) as a function 
of (k x /n,k y /n) for the UICtItIt 1 state. 



FIG. 6: ((a) Contour plot of the spinon dispersion E+{k) 
as a function of (k x /n,k y /n) for the Z2At°t+t+ state, (b) 
Lower edge of the two-spinon spectrum E2 S (k) as a function 
of (k x /n,k y /n) for the Z2At°t°t+ state. 



liquid), which is given by 

Ui, i+ x = X t1 , 

Ui,i+y = XT 2 , (28) 

1,2,3 n 

The spinon dispersion is 

E+(k) = 2xy / cos 2 k x + cos 2 k y (29) 

The spinon dispersion E + (k) and the lower edge of 
the spin-1 spectrum E2 S (k) are plotted as functions of 
(kx/n, ky/n) in Fig. 5. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -OS 0.0 5 1.0 



(a) k < / " (b) 



The second one is the Z2At°t+t+ spin liquid state in 
Ea. (|15|) . Energy dispersion of the Z2At°t ( It1 state is 
given by 

E + (k) = + 

£ i = X ( cos k x + cos k y ) + X cos (k x + k y ) — a\ 
e 2 = V {cosk x — cos k y ) (31) 

The spinon dispersion E + (k) and the lower edge of 
the spin-1 spectrum E2 S {k) are plotted as functions of 
(k x /w,k y /Tr) in Fig. 7. 




(a) (b) 



FIG. 5: a) Contour plot of the spinon dispersion E+(k) as a 
function of (k x /n, k y /n) for the 7r-flux state, (b) Lower edge of 
the two-spinon spectrum E'2 3 (k) as a function of {k x /n, ky/n) 
for the 7r-flux state. 

Now we consider two other Z2 spin liquids. The first 
one is the Z2Ar r"r^ spin liquid state in Eq. ijTlJl . The 
energy dispersion of the Z2At°t+t!1 state is given by 

E + (k) = 2^/ej + el 

e i = X ( cos k x + cos k y ) + X cos (k x + k y ) — a\ 
£ 2 = V (cos foe + cos k y ) (30) 

The spinon dispersion E + (k) and the lower edge of 
the spin-1 spectrum E 2s {k) are plotted as functions of 
(kx/ir, ky/n) in Fig. 6. 



FIG. 7: (a) Contour plot of the spinon dispersion E+(k) 
as a function of {k x /n, k y /n:) for the Z2At x t\t\ state, (b) 
Lower edge of the two-spinon spectrum E2 3 (k) as a function 
of (k x /n, ky/n) for the Z2At 1 t]_t!1 state. 

If the mean-field state is stable against the gauge fluc- 
tuations, we expect that the mean- field spin-1 spectrum 
E2 S should qualitatively agree with the real spin-1 spec- 
trum. 



VI. SPIN SPECTRAL FUNCTION 

The quantity of interest for comparison with experi- 
ments is the spin spectral function in (q, oj) space S(q, w). 
S(q, u>) can be calculated from the Fourier transfor- 
mation of two spin correlation function in real space 
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(Si(t) ■ Sj(0)}. Due to the spin-rotation invariance, we 
obtain that 

(Si ■ Sj) = 3(Sf ■ 5?) = 3(5f • Sf) = 3(Sf ■ 5J> 

The spin spectral function at zero temperature is given 
by 



S(q,u)=3S"(q,<j) 



(32) 



where 



S zz (q, u)=J2 27r l < A I SZ (l) \G) \ 2 S(E G +u;- E x ), (33) 



S Z {q) 



E 

j 



2 ^ / < ^k+q 
k 



the summation is over all the eigenstates |A) and N is 
the number of lattice sites. 

Numerical calculation indicates that the neutron scat- 
tering spectrum predicted by the dynamics spin correla- 
tion contains continuous modes besides the main sharp 
peak. The Fig. 8-11 show the spin spectral function 
S zz {q, u) in the q-space. The point T is (0, 0), M is (0, tt), 
and X is (7r,7r). Darker shade represents a larger spec- 
tral function. These results can be compared with the 
inelastic neutron scattering experiments. The scan from 
r to X corresponds to the scan from k — to k — 2n/b 
in Fig. 2a of Ref. [3. 



-e 
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FIG. 9: The spin spectral function S(q, uj) of the U1Bt x t^t\ 
state in Eg. (1221 . The solid lines mark the lower edge of the 
spectrum. 




FIG. 10: The spin spectral function S^,^) of the IT CtJt-T 1 
state in Eg. 1231 . The solid lines mark the lower edge of the 
spectrum. 



FIG. 8: The spin spectral function S(q,u) of the SU°°(2) 
spin liquids in Eo. 1240 . 



VII. COMPARISON WITH EXPERIMENTS 
AND PREVIOUS THEORETICAL RESULTS 

Our J- J' model is designed to describe the spin fluc- 
tuations in the CS2CUCI4 sample. The spin dynamics in 



CS2CUCI4 at low t emp eratures were measured by neu- 
tron scattering. To the first order approximation, 
the measured spin spectral function is similar to that 
of decoupled ID spin chain (or the SU°°(2) spin liquid) 
(see Fig. 8). The decoupled ID spin chain has gap- 
less spin-1 excitations along the lines k x + k v = ir and 
k x + k y = 2ir. If we examine the measured spectral func- 
tion more closely to determine the low energy bound- 
ary of spin-1 excitation £i(k x ,k y ), we find the gapless 
spin-1 excitations actually only appear at certain isolated 
points, implying that the spin state is 2-dimensional. The 
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FIG. 11: The spin spectral function S(q,u>) of the n -flux 
phase (SU2B state) in Eg.i25t 

experiment 0,0] strongly suggests that there is a gap- 
less spin- 1 excitation at (k x ,k y ) = (ir/2+e, ir/2+e) where 
e ~ 0.04. However, it is unclear whether the spin-1 ex- 
citation near {k Xl k y ) = (tt,tt) is gapless or has a small 
gap. The experiment also determines the high energy 
boundary S2(k x ,k y ) of spin-1 excitations. It was found 
that £2(k x ,k y ) has maxima near (k x ,k y ) = ±(7r/2,7r/2) 
(the location is not known precisely). 

Using this information, we can determine which spin 
liquids classified in this paper cannot describe the spin 
liquid state in CS2CUCI4. First the 7r-flux phase has no 
low energy spin-1 excitations near (7r/2, 7r/2). Thus it 
cannot describe the spin liquid state in CS2CUCI4. The 
U1Bt 1 t^_t\ also does not work since the spin-1 excita- 
tions open a small gap near (jr/2, 7r/2). 

It appears that the UICt+t^t 1 spin liquid and some 
of its close relatives such as Z2At°t+t+, Z2At 1 t\t\, 
Z2At 1 t 3 _t s L and Z2At°t^t^_ spin liquids are consistent 
with observed spin spectral function. Let us list the 
universal characters of the above five spin liquids. The 
UICt+t^t 1 and Z2At q t\t\ spin liquids only have gap- 
less spin-1 excitations at (0,0), (7r,7r), ±(7r/2 + e, 7r/2 + e) 
and ±(tt/2-£, ir/2-e). The Z2At x t\t\ spin liquid only 
has gapless spin-1 excitations at (0, 0), ±(n — e', n — e'), 
±(7r/2 + e,7r/2 + e) and ±(tt/2 - e,7r/2 - e). The 
ZIAt rt-T^L and Z2At°t}_t'^ spin liquids only have gap- 
less spin-1 excitations in small patches near (0, 0), (77, 7r), 
and ±(tt/2, 7t/2). By measuring the locations of gapless 
spin-1 excitations and their spectral functions, we hope 
further experiments can determine which spin liquid ac- 
tually describes the sample. 

However, to compare the spin spectral function quan- 
titatively with experimental measured results. |l8tll9l| we 
need to include fluctuations around the mean-field state. 
The interaction due to fluctuations enhances the spin 
fluctuations at (ir/2,ir/2). The strong enough interac- 



tion can drive the transition from the spin liquid state to 
the spiral state. 

To include the interaction enhancement of the 
(7r/2,7r/2) spin fluctuations, we consider the effects of 
the J-term (the diagonal coupling term): ^\ JSj-Sj +x+y . 
Within the random phase approximation (RPA), the spin 
correlation function 7r(k, ui) can be obtained from the 
mean- field spin correlation function 7r (k, ui): 



1 — J cos(k x + /c a )7To(k, ui) 

The imaginary part of 7r(k, ui) gives us the RPA spin 
spectral function. To be definite, we start with the mean- 
field UICt+tS. t 1 spin liquid state and choose 7To(k, ui) 
to be the mean-field spin correlation function of the 
UICt+t^t 1 state. The experiments 0,0] indicate that 
the system is close to the spiral phase. So here we choose 
the value of J to make the system close to the spiral 
phase instability. The resulting RPA spin spectral func- 
tion is plotted in Fig. 1121 and Fig. ^| The spin spectral 
functions along the lines A, C, and D in Fig. ^| were 
measured by experiments. |l8j The simple RPA results 
are quite similar to the measured results (see Fig. 3 in 
Ref. 0). 
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FIG. 12: The spin spectral function S(q, ui) of the UICtItIt 1 
state obtained by RPA. The short solid lines mark the lower 
edge of the spectrum. 

We would like to remark that the low energy effec- 
tive theory for the UICt+t^t 1 state contain gapless 
fermions interacting with U(X) gauge field. Depending 
on that strength of the U(l) gauge fluctuations, the 
U(l) gauge interaction may or may not destabilize the 
UICt+t^t 1 state. It is shown that if the spinons has 
a very anisotropic dispersion, the instanton effect of the 
2+1 U(\) gauge theory is irrelevant. pj In this case, the 
UICt^t^t 1 state is an algebraic spin liquid. |30l| which 
contains gapless spin-1 excitations, but non of the gapless 
excitations are described by well defined quasiparticles. 
If the spinons has a more isotropic dispersion, the instan- 
ton effect is relevant, which lead to a confinement of the 
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FIG. 13: The RPA spin spectral function S(q, ui) along the 
lines A, C, D, G, and J in Fig. [Hfor the UICtItIt 1 state. 

U(l) gauge field. In this case, the U1Ct+tH t 1 state may 
become one of the stable Zi spin liquid states or other 
symmetry breaking states (such as the spiral state) below 
a transition temperature T c . 

Even when the UICt+t^t 1 state is unstable, it is still 
useful for understanding the system at finite tempera- 
tures. The UICt+t^t 1 state may control the dynamical 
properties of the CS2CUCI4 samples above the T c , pro- 
vided that T c is much smaller than the band width of 
the spin fluctuations. 



VIII. CONCLUSION 

In this paper, we use the projective symmetry group 
(PSG) to characterize quantum orders in symmetric spin 
liquids on 2D triangular lattice. We classify these sym- 
metric spin liquids through PSG, and find that there are 
63 kinds of Z 2 spin liquids, 30 kinds of U(l) spin liquids 
and 2 kinds of SU(2) spin liquids. The mean-field phase 
diagram Fig. 2. for the nearest neighbor spin coupling 
system is calculated. One-spinon and two-spinon exci- 
tation spectrum for some Z2 and U(l) spin liquids are 
inspected. We show that the gapless spinon-1 excitation 
spectrum can be used to physically measure the quan- 
tum order. By examining the measured spin spectral 
function in CS2CUCI4 sample. [Tsj we find that it is pos- 
sible that the spin liquid state in CS2CUCI4 is described 
by the U1C r^r^T 1 spin liquid or one of its relatives. 

We would like to point out that the RPA spin spec- 
tral function for the U1Ct+ t° t 1 spin liquid is very close 
to the spin spectral function 25] obtained from the slave- 
fermion approach. [2^ In the experiments. [Tsl the 
spin liquid was observed only at a low but finite temper- 
ature. At that temperature, it is hard to distinguish to 
the two proposed states. However, if we find a spin liquid 
state at zero temperature, then we can tell which of the 
two proposed spin liquids can describe the experiments. 
This is because at zero temperature, the U1Ct+t2_ t 1 spin 
liquid (and its relatives) and the spin liquid from the 



slave-fermion approach are qualitatively different. The 
UICt+t^t 1 spin liquid (and its relatives) contains gap- 
less spin-1 excitations while the spin liquid from the slave- 
fermion approach is fully gapped. 

XGW is supported in part by NSF Grant No. DMR- 
01-23156 and by NSF-MRSEC Grant No. DMR-02- 
13282. 



APPENDIX A: CLASSIFICATION OF Z 2 SPIN 
LIQUIDS 

In this appendix we consider the connected ansatz only, 
which is interesting and cover a wide range of spin liquids 
in fact. For any two sites i and j, a connected ansatz will 
offer a path {i\i2 • • - in), which satisfy that Ui k ^ k+1 ^= 0, 
k = 1,2,..., 71, i\ = i and i n = j. Due to the transla- 
tion symmetry of the ansatz, we can choose a gauge in 
which all the loop operators of the ansatz are transla- 
tion invariant. We will call such a gauge uniform gauge. 
Hereafter we will work in uniform gauge. According to 
the definition of PSG and uniform gauge, 

G x (i)P TxC Gl(i) = G x {i)P c Gi(i) = Pc, 
G x {i)P c = P C G x (i). 

For Z2 spin liquids, different loop operators basing at the 
same base point do not commute. We see that translation 
invariance of Pc in the uniform gauge requires that 

G x (i),G v (i)eG = {±T }. 

Gauge transformations W(i) = (— )^ r° (f(i) = ±1) do 
not change the translation invariant property of the loop 
operators. We can always choose a gauge which satisfy 
that : 

G y (i)=r°. (Al) 

Gauge transformation W(i) = W(i x ) does not change 
the condition G y (i) — r , hence we can gauge fix that 

G x (i x ,i y = 0) = T°. (A2) 

From the relations (|13|) . it is easy to see that 

G x (i) G y (i - x) G- 1 (i - y) Gy 1 (i) = ±r°. (A3) 

Substituting Eqs. (|A1|) and l|A2(l into the above, one sees 
that 

G x (i)G x 1 (i-y)=±r°. (A4) 

Therefore there are two classes of PSG belong to trans- 
lation symmetry, 

G x (i) = G v (i) = t° (A5) 

and 

G x ( l )^(-)^r°,G y ( l )^r . (A6) 
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We will label the former Z2A and the later Z2B. 

Then we consider the spin parity symmetry T* which 
satisfies that 

T * TxT *-i T -i = T*T y T*~ 1 Ty 1 = T* 2 = ±1. (A7) 

In a similar way, we have the following relations 

Gt^GMGt-T*)- 1 ^^)- 1 e Q, 
Gt*T* G y Ty{GT'T*)^ 1 {GyTy)^ 1 e G, (A8) 
G T *T*G T *T* e G, 

then 

G T * (i) G x (i) G T » (t - i) G- 1 (i) = r^r , 

G T . (i) G B (i) G T i (t - y) G- 1 (i) - %t r°, (A9) 

Gt* - ±T°, 



Gt* (i) = r)i%V y v t9T*, 



(A10) 



where = ±1, rj yt = ±1 and g\, — ±r°. Two gauge 
inequivalent choices of t/r* are gx* = t° and 57-* = «t 3 . 

Next we add two types parity transformations 
Pxy [ixjiy) — (iyjix) and P X y [ixiiy) = ( — iyi ~ix)i which 
satisfy that 



J-x^xyJ-y r X y 


= 1, 


± y r xy ± x r xy 


= 1, 


T P -T P - 

± x 1 xy- 1 y 1 xy 


= 1, 


T P -T P - 

J -x 1 xy 1 - y 1 xy 


= 1, 


rp* p rp* — l p— 1 

r xy 1 r xy 


= 1, 


rp* p rp* — 1 p 
^ r xyl- r xy 


= 1, 



(All) 



From the above equations, we find that 

(G x T x )(Gp xy P xy )(GyT y y 1 (Gp xy P xy y 1 e G 

(GyTy)(Gp xy P xy )(G x T x ) 1 (G Pxy P xy ) 1 e G 

{G X T x ){Gp :l;y P X y){GyTy){Gp :EV P X yy 1 £ Q 

{GyT y ){Gp xv P X y){G x T x ){G Pxv P X yy 1 £ G 

GT'T\Gp xy P X y){G T ,T*Y\G Pxy P X y)- 1 e G 

GT*T*(Gp xy P X y)(GT»T*)~ 1 (Gp xy P X y)~ 1 e G 



(A12) 



and 



where rj xpxy = ±1... 
For Z2A type PSG's, 

^ ' Pxy (0 VxpxyVypxy9P xy 7 

Gp x y(i) = VxpxyVypxygPxy- (A14) 

For Z2B type PSG's, 

Gp xy («) = (-) % *" rtxpxyVypxydP.y, 

Gp x§ (l) = (-f^o VxpxyVypxy9P X y- ( A1 5) 

We also have that 

^xj/ = Pxy = PxyPxyP xy P X y = 1 (A16) 

which leads to 



G PsC y^)G Pxy {P X y{i)) = ±r°, 
Gp. s (i)Gp xS (P x ,(i)) = ±r°, 



and 



±r°, 



(VxpxyVypxy) V 9p xy 
(VxpxyVypxy) V 9p xy = ^ szT J 

The above implies that 

Vxpxy — Vypxy — Vpxy? 
Vxpxy — Vypxy = Vpxy? 

g Pxy = ±t\ 9 % v = ±t°. 
By the similar way, one sees that 

Vxt = Vyt = T] t . 



(A17) 



Gauge transformation W(i) = {—) %!C will change rj pxy to 
—Vpxy, therefore we can chose a gauge in which r\ pxy = 1. 

In conclusion, we have the following 60 kinds of Z2A 
algebraic PSG's and 60 kinds of Z2B algebraic PSG's on 
our triangular lattices: 

Z2A 9p xv {gp X y )r, pX y (gr* )n t 



G X (i) = Gy(i)=T , 

Gt* (i) = vl9T*,Vt = ±1, 
Gp xy {i) = gp xy , 



(A18) 



Gx (*) G Pxy (i - x) G" 1 (P xy (*)) G P l (i) 


— ^]xpxy7~ 


Gp„ s (i) 


2 1 -i 

— VpxygPxyiVpxy — ^tJ-j 


G„ (*) G PxH (i - y) G" 1 (P xy (*)) G P l (i) 


= T }ypxy T 


and 




Gx (i) G PxS (i - i) G" 1 (i) - 4) G Px x e (i) 


T]xp%yT 


Z2B gp* v {gp*y)n P , 




G v (*) G Pxe (t - y) G- 1 (i) - y) G^ s (i) 


= Vypxy 7 ' 


Gy{i) = 


r°,G x (i) = (-)'» r°, 


Gt* WG Pxy (i)G T l(Px V (i))G P l v (i) 


~ VtpxyT 


Gt* (i) = 


vl9T',Vt = ±1, 


Gt* WG Pxe ( J )G T , 1 (P^(*))G Px 1 e ( J ) 


— T]tpxy7~ 


Gp,y{i) = 






(A13) 


G Pxi {i) - 





(A19) 
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where g T * , gp xy , and g Pxy satisfy that 



Ifrp* 


= ±T° 




= ±T° 


z 


= ±T° 


9T-gp xy 9T^9pl y 


= ±T° 


9T*gp xy gT*9p* y 


= ±T° 


HI' ,9P .//;■'.'/;-;. 


= ±T° 



All the gauge inequivalent gr* , gp xy , and gp xS are given 
by the following table 



<7t* 


Qp 

■y i xy 


gp x ~ 


Qp Qp - 


T° 


T° 




T° 


T° 


T° 


ir 3 


IT 3 


T° 


IT 3 


T° 


IT 3 


T° 


ir 3 


ir 3 


T° 


T° 


IT 3 


ir 1 


IT 2 


7T 3 


T° 


T° 


T° 


ir 3 


T° 


ir 3 


■ 3 

IT 


7T 3 


T° 


ir 1 


ir 1 


IT 3 


ir 3 


T° 


IT 3 


ir 3 


IT 3 


ir 3 


T° 


ir 3 


■ 3 
IT 


ir 1 


IT 2 


ir 3 


ir 1 


T° 


IT 1 


ir 3 


ir 1 


ir 3 


T 2 


ir 3 


ir 1 


ir 1 


T° 


ir 3 


ir 1 


IT 2 


IT 3 



After finding all the PSGs, next we need to find the 
ansatz that is invariant under those PSGs. This way we 
classify (mean-field) Z 2 spin liquids through PSG's. Let 
us first consider the translation symmetry. We have two 
classes of spin liquids, 

Z2A spin liquids 

Gy(i) = T Q ,G X (l)=T°, 

(A21) 

and 

Z2B spin liquids 

G y (i) = T°,G x ( l ) = (-y«T°, 

Ui, i+m = (-) m ^u m . (A22) 
Considering the spin parity symmetry T*, we have 

G T *T* (u M+m ) = u iii+m , (A23) 

which reads 

- n^g T *u m g^l = u m (A24) 



Therefore, for different {gr*)^, we have the following 
table: 





u m = 




u m — 0,if m — even 


4 


u m = u m T X + u 2 m T 2 


T 3 


u rn = u m rl + u m r2 i ^ m = even 

u rn = u m r° + u 3 n T 3 , if 772 = odd 



Now we consider the 180° rotation symmetry, which 
can be constructed through combining P xy and P xy , 

Gp^PxyGp^-Pxyiui^i+m) = Ui^ +m (A25) 
the above leads to 

vZy 9P xy 9P xy u m 9p] y 9~ P ] y = ul , (A26) 
for the Z2A spin liquids, and 

R m * m " vZ9P xy 9P xV u m g^g- p l = u m . (A27) 

for the Z2B spin liquids. 

By writing all the ansatz, we find that although there 
are 120 kinds Z 2 algebraic PSG's, only 23 of them lead 
to Z2A spin liquids and 40 of them lead to Z2B spin 
liquids. Other algebraic PSG's lead to vanishing ansatz, 
17(1) or S77(2) spin liquids. 

APPENDIX B: CLASSIFICATION OF U(l) SPIN 
LIQUIDS 

For U(l) spin liquids, we can choose a gauge so that 
Uij take the form 

Uij = ip t3 e l6 ^ 3 = u%t° + u%t 3 (Bl) 

We will call this gauge canonical gauge. In the canonical 
gauge, IGG has a form of Q = je ier '\ 6 G [0, 2tt)} , where 
8i is real for each i. Due to the translation symmetry of 
the ansatz, the loop operator Pq have a form 

Pt x c=(tT P c(tT (B2) 

where rn = 0, 1. The gauge transformation G X:V associ- 
ated with the translation take the following form 

G x (i) = Hr 1 )"*^^ 3 (7T 1 )"-* 

G v (i) = (-zt 1 )" 1 e 16 ^ (ir 1 )™ 1 "* (B3) 

in the canonical gauge. We note that a gauge transfor- 
mations keep to have the form in the canonical gauge 
must have one of the following two forms 

Wi = e i9 ^ 3 (B4) 
Wi = e 10 ^ 3 (ir 1 ) (B5) 
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For spin liquids with connected Uy, G XtV must take one 
of the above two forms in the canonical gauge. Thus rij 
can only be one of the following four choices: n, = 0, rij = 

(i - (-) l )/2, m = (i- (-HA and m = (1 - (-H/2. 

In these four cases, G^y take one of the above two forms 
and Uij can be connected. 

Let us consider those cases in turn. We will work in 
the canonical gauge. 

When rii = 0, G x ^ y have a form 

G x (i) = e^ Wr3 

G y (i) = e i6 y^ 3 (B6) 

by (|B4|I we can gauge fix them up to 6 X (i x , 0) = and 
y {i) = 0. The relation (|T3*)) requires that 

X (i) + y (i - x) - 0* (i - 0) - y (i) = ip (B7) 

Therefore G K , y have a form 



Using gauge transformation W% — l »v/ 2 ) T j the above 
change to 



G x (<) = e^ T (Jr 1 ) 

= e"-* 8 (B16) 

The ansatz with translation symmetry has a form 



When n -, 



(l-(-)%)/2, 



i(-) !x 0„ 



(B17) 



then 



G x (i) 
G„(i) 

G x (i) 
G„(<) 



s i((-)'« V +9 x )r 3 



G x (i) = ei^v+s-)^ 3 
G„(i) - e^ S 



(B8) 



Using gauge transformation Wj = e l "~ 
above change to 



(B18) 



(B19) 

"V/2)r 3 the 



The ansatz with translation symmetry has a form 

u iti+m = ip m e^- m " i ^ + ^ r3 (B9) 
When m - (1 - (-) 4 )/2, 

G x (i) = e^^iir 1 ) 



G y {i) = e^W^^r 1 ) 



(BIO) 



then 



G x {i) = e i6 * T 

G y (i) = e M » r '(ir 1 ) (B20) 
The ansatz with translation symmetry has a form 

"M+m = ip m e l( - yv ^ r3 (B21) 

Now we add parity symmetry P xy . Considering the 
flux in a plaquette, the parity symmetry requires tp = 
or 7r in (IB9II. From the relations 



G x {i) = e *((-)'»*'+fl-)'- s (i T 1 ) 




G, (i) G PxB (i - sc) G" 1 (P XJ/ (<)) G^ s (i) 


7l9 T 3 

_ — £ xpxy 1 


G tf (i) = e^Cir 1 ) 


(Bll) 


G, (i) G Pxs (i - y) G- 1 (P xy (*)) G p i y (i) 


if) T 3 
— p La ypxy ' 


Using gauge transformation VU; = e*^ - ' " 


"V/^r^ the 


G, (*) G PiS (i - i) G; 1 (P xy (i) - i) g£ s (i) 


if) -t 3 
— ^ LU xpxy ' 


above change to 




G, (i) G PiS (t - y) G- 1 (P^ (i) - y) G?\ (i) 


if) -t 3 
— g LU ypxy I 


G x (i) = e i9 - T '(ir 1 ) 




G Pxv (i) G Piy (P xy (<)) 


if) t 3 

— giUpxyl 




(B12) 


Gp. s «G PiS (P xy (i)) 


fJ-QpxyT 


The ansatz with translation symmetry has 


a form 


Gp, y (i) G PxS (P XJ/ (i)) G^ (P xg (*)) G^_ (i) 


= e i6 v r 3 



When n< = (1 - (-) 2 *)/2 



»(-)Vn 



(B13) 



then 



G x (i) = e^-W^fir 1 ) 



(B14) 



(B15) 



(B22) 

we will find only rii = and n.j = (1 — (— ) l )/2 lead to the 
existence of Gp^ (i). Then we can classify the U(l) spin 
liquids with parity symmetry into three classes according 
to their translation property. 
UlA spin liquid: 



G x (i) 

Gyii) 



e x 

e i0yT 3 

ip m e 1 ' 



(B23) 



14 



UIB spin liquid: 

G x {i) 

Gy{l) 

U1C spin liquid: 



_ e i(i'y7T + ex)l- 3 

= iH 



(B24) 



G s (i) = e^dr 1 ), 

,»(-)> m r 3 



(B25) 



ip m e 



Next, we should find all the gauge inequivalent PSG's 
Gp xy (i) and Gp x9 (?) for the above cases respectively. In 
fact, as proved in Ref.Q, it exists the following relations 
between UlA and UIB PSG's: 



Gt* {i) = G T * (<) 



(B26) 



where Gp XH! p x9 are UlA PSG's. Hence we can consider 
the UlA PSG's only, and using the relations, we can 
construct the UIB PSG's. 

1) UlA spin liquid, G B (») = e t6 * T \ G v (i) 
Ea. (|B22)l reads 



and 



Gp xh (i) 
or Gp xh (i) 

Gp xS (i) 
or Gp xS (?) 



y ) 'Ppxy ~l~yp X y 



4>p 



(B27) 



— z>^((^aj ""H^y 'j^Ppxy "H^pxi/ ) 7 



gt((*aj iy)<£'pxy~\~4 ) pxy)'T ( ^^-l^j 



Gauge transformation W(i) = e *(. % *T*v'i<PpxyT 3 /2 t ra ns- 
fers Ell and dB28l) into 



Gp xh (i) 



— p^'Ppxy 



T " or G Pxh (i) = e *«W 3 (zr 1 ) , (B29) 



and do not change the forms of G Xjyi p xS . From the last 
equation of (|B22(I . we obtain the following four cases, 
a) 

Gp xh W = e < W\ 



Gp xS (») 



~^~^"y ) < ftpxy~\~4*pxy') T 



Gauge transformation = e i ( t z+ 4 »)</>p*s/ 2r3 trans- 

fers the above into 



Gp xh (i) 
Gp xS (i) 



(B30) 



and does not change the forms of G x y . 
b) 



Gp xy (i) 

Gp xS W = '/;.,V"' ; '"(^' )•'/;.., tl (B31) 



c) 



Gp xh (i) 
Gp xS (i) 



Vpxy e ^ PSCV ■• Vpxy = ±1 



(B32) 



d) 



GPxy (t) 
GPxy (i) 



£>^{ (^flS ^y ) ^Ppxy ~~\~*Ppxy ) T 



Gauge transformation W(i) = e l s/)¥w/ 2 ' r3 trans- 
fers the above into 



GP x y (i) 

Gp xS (i) 



(B33) 



and does not change the forms of G x v . 

Now we add the last symmetry spin parity T*. For the 
same reason, Gp* (i) can take one of the two forms l|B4|) 
and <|Bf>|) . The relations (|A12|I require that 



G T - (*)<?„ (i)G^(i-x) G- X (i) 
G T , (i)G y (fiG^ii-fiG- 1 ® 
Gt* {i)G Pxy (i)G T l(P xy {i))G- p ] y (i) 
Gt* (»)Gp. s (i)G^(P xS (i))G^ 9 (i) 

G|. (i) 



Hence we have that all the UlA type PSG's: 



p^"tpxy' 



pL"tpxy> 



(B34) 



GPxy (*; 


I — pi^pxyT 




GPxy (»; 






Gp. (i) 




{it 1 


GPxy 


]\ = p^^pxyT 




GPxy (i 


= ^e'W (V) 




G T . (i 


5) = ijje*** 3 e^'l 




GPxy 


i) = e^i-^ 3 (rr 1 ) 





(B35) 



(B36) 



GPxy (*) 

G T , (*) 



i _ i<t>p Xy T 



(B37) 
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Gp xy (<) = e*™** (zr 1 ) 

Gp XH (*) = e^ 9 " 3 (ir 1 ) (B38) 
Gr.(») = ^e^k^-i^je^^ir 1 ) 

We can label these spin liquids as 
UlAgp^igp^^igT*)^, where g Pxy , Pxy = r ' 1 
corresponding to e i * px «- pa: s' r and e* I! '' , " : » T (* t1 )- Now 
we consider the form of ansatz that in invariant under 
the above PSG's. The translation symmetry requires 
that 



(B39) 



The 180° rotation symmetry requires that for 



Other PSG's lead to vanishing ansatz or SU(2) spin 
liquids and can be dropped. 

2) U1B spin liquid, G x (i) = (-)*" e l9 * T \ G v {i) = 
e iB y T _ p rom relations l|B26(l . one can obtain all the U1B 
PSG's easily, 



Gp„,(i) = 
Gr.(i) : 

Gp, y (*) 
Gp, 9 (*) 
G T - (<) 



(-r ' 



(B47) 



(B48) 



and for gp xy {gp xt ) 



r o T i 10 



Vpxy ' 



(B40) 



Gp.„(<) = (-rH^e^ 3 (B49) 



W m — 1 lpxy T u -rn T — Vpxy U m 



(B41) 



The spin parity symmetry T* requires that for 


Gp. y (<) 


= (-r 


x«y gi^pi 


^ (zr 1 ) , 




Gp, e (<) 


= (-)' 


xiy e i(ppoL 






Gt- (*) 


= ^ 


\Vt=- 




U m = - Um (B42) 


We can 


label 


these 


spin 



and for (<?t*)? 



(B43) 



We find that the following 6 sets of ansatz that give 
rise to U(l) symmetric spin liquids: 



(B50) 



Now we consider the 



UlBgp^gp^^gT*)^. 
form of ansatz that in invariant under the above PSG's. 
The translation symmetry requires that 



(-)' 



(B51) 

The 180° rotation symmetry requires that for 



u- -.L. - u° r° + u 3 t 3 



,0,3 



0, if to = even 



(B44) 



Here we have used the notation U lA[ab, cd]e to represent 
the collection UlAabe and UlAcde. 



J71 J 4[r r?,r 1 rI]rI 



UlA^T^T 1 ^]^ 



vLt 3 



0, if m — odd 



(B45) 



(B46) 



u m = (-) 

and for g Piy (gp x 



m x m v f 



(B52) 



r o T i r i T o 



u m = -(-) m * mv r 1p n X yUm (B53) 
The spin parity symmetry T* requires that for 

(9T*) nt =T°_ 



and for (gx 



* )Vt — T T) t 



(B54) 



(B55) 



The ansatz gives rise to U(l) x U(l) spin liquids since 
Uij only connect points within two different sublattices. 



We find that the following 10 sets of ansatz that give 
rise to U(X) symmetric spin liquids: 
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[/1B[t tA,tV]t° 

u m 3 = 0, if m = even 
UIB[t tI,t 1 tX]t\ 

Ui,i+m = (-) <r 



(B56) 



0, if m x = odd and m y = oc?<i (B57) 



u- -.l. - f-^ 1 - u 3 r 3 

v? m = 0, if m x = even and m y = even (B58) 

UIB^t^t^t 1 

0, if = odd or m y — odd (B59) 



H,i+m — (~ ) " 7 ' 

u 3 



only connect points within four different sublattices. 
I71B[T°Ti,r 1 r^]Ti 

-M +m = (-)*"«*• 

u hi = 0; if 771 = even 

uf n = 0, if = even or m 9 = even (B60) 



Other PSG's lead to vanishing ansatz or SU(2) spin 
liquids and can be dropped. 

3) U1C spin liquid, G x (i) = e* 6 *^ (it 1 ) and G y (i) = 
e ie y T \iT l ). 

G Pa!y (i) = e'W^e'W^ (ir 1 ) . 
For the former four equations in (|B22|> . we have that 

@pxy(i) ( ) Ppxy ~t~ 4*pxy: 

ix^ ~t~ ( ) (Ppxy ~t~ 4*pxy 

and the last three equations in (IB22I) require that 

r?n — n l n ix P i ( ? 'xpx ! , 7r / 4 +^y) 7 " 3 

^Pxy \ L ) — 'Ipxy'lxpxye ' ! 

g«((-)Vpxy+0px!,)T 3 (j-,- 1 ^ (B61) 

and 

I'J — 'Ipxy'lxpxye ' i 

where r)p X y,xpxy = ±t,Vpx y ,xpxy = ±1. Now we add the 
last symmetry spin parity T*. In a similar way, through 
the relations (|A12|) . we have that 

G T * (i) = ryje*- 3 , 

^x e <(-)V t +^ t )- 3 (fri) (B63) 



where ?7 t = ±1,77^4 = ±1. 

It is noted that gauge transformation Wi = e 1 ^ ^ T 
and Wi — njp^y do not change G XtV , and can be used to 
simplify the forms of Gp x : p xgt T*- We conclude all the 
U1C type PSG's: 



G P 

^xy 


(i) 


Vxpxy* 


G P - 

x y 


(0 


Vpxy^l 




(i) 


= Vie 1 *' 



Gp, y ii 
Gp* y (i 
G T * (i 



Uij 


Gp 

*xy 


(i) 


= e <^xp 




G P,y 


(i) 


Vpxy^l 




Gt* 


(0 


= 4^ 



G P (i 

r xy \ 

G P - (i 

ixy \ 

G T * (i 



g l (')xpx !/ 7r / 4 +^px!/)-r 3 

^ e «(^xpx y -/4+0Px 9 )r 3 (B64) 
pxy^ 



I 1 n l * e ilptT (ir 1 ) 

|rji = — li 'lxpxy c V I 



= ri 1 * p*(''xpx H 7r / 4 +^PX !/ )T 3 

Ixpxy 

= n l vxy e l ^v^l^^y (it 1 ) (B65) 



ipxy 



(iT 1 ) 



f*e<^/ 4+ ^y (B66) 



e ^xpxW4+0px y )r 3 ^ 

n pxy e l ^y* /4+ ^y (it 1 ) (B67) 



i i(j>tT I 



(iT 1 ) 



We can label these spin liquids as 
U^C(gp scy ) ri ^ y (gp xi ) llpsc9 (gT*)- Now we consider 
the form of ansatz that in invariant under the above 
PSG's. The translation symmetry requires that 



7*i,i-f-r. 



u° m T° + (-yui 



(B68) 



The 180° rotation symmetry requires that for 

(;,.,, [i)G,:Sr, :r n) 

u m = ~n pX yU m: u m = (—TJpxy) v m ] (B69) 

for GPxy (i)G Pxy (Pxy(i)) = VpxyVxpxy (iT 1 ) 



u° = -n m n my u° 

in Ipxy Ixpxy "'m' 

u in = ^(^Vpxy) Vxpxy u m- 



(B70) 



The spin parity symmetry T* requires that for 



G T * (i) = n\e^ 

_ ( _\ m 



(B71) 



and for G T * (i) = *7jj« v e^ T (ir 1 ^ 



v° = -n m * v° v 3 = 77 mx 7y 3 CB721 

'Ixpxy "mi ' Ixpxy y 1 -" ^ ) 
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We find that the following 14 sets of ansatz that give 
rise to U(l) symmetric spin liquids: 



Ui, i+m = U° n T° + (-)* U 3 m T 3 

u' = 0, if m = even; 



UlC[rlrl,rlriy 



i-Yulr 3 
0, if m — odd 



(B73) 



(B74) 



(B75) 



The ansatz gives rise to £7(1) x U(l) spin liquids. 

u i,i+m — ( — ) u mT i 
^i,i+m ( ) ^m^~ 

u m = 0, if = odd or m y — odd; (B76) 
The ansatz gives rise to (t/(l)) 4 spin liquids. 



u m — 0, if m x — even or m y = odd 
v? m — 0, if m x = odd; 

0^1 „-1^0l„-l 



(B77) 



U1C[t!LtI,tIt*L]t 



u° m r° + (-)< u 3 m r 3 



u® n = 0, if = even 

v? m = 0, if m x = odd or m y = even; (B78) 

Other PSG's lead to vanishing ansatz or SU{2) spin 
liquids and can be dropped. 



And we can gauge fix the above form as 

G x (i) = 9x 
Gy{i) = Vy(i)9y 



(C4) 



where r) y (i x = 0, i y ) = 1. Now the relation l|13fl reads 

r, y (i-x)r, y (i)€SU(2) (C5) 

We find that there are only two different PSG's for trans- 
lation symmetric ansatz 

G x (i) = g x ,G y (i)=g y (C6) 
and G x (i) = g x , G v (i) = (-) 1 * g y (C7) 

The two PSG's lead to the following two translation sym- 
metric ansatz 



and u i>i+m = (-) 



(C8) 
(C9) 



Similar to Zi case, adding more symmetries, we have all 
the £{7(2) PSG's belonged to triangular lattices. 

SU2Arj pX yrj t 



G x (i 

Gy{l 

G T *{i 
Gp xy {i 
G Pa M 



9x, 
9y, 

Vt9T* , 

Vpxy9 Pxy 1 



(C10) 



and 

SU2Bri pxy rit 



APPENDIX C: CLASSIFICATION OF SU(2) SPIN 
LIQUIDS 

We assume that for a SU(2) spin liquid, one can always 
choose a gauge so that Uij has a form 



u i:i = u^t 



(CI) 



Hereafter we call this gauge canonical gauge. In the 
canonical gauge, IGG has a form Q = SU(2). The gauge 
transformations that keep Uy to have the form in the 
canonical gauge are given by 



Wi = rj(i)g 



(C2) 



where ?y(i) = ±1 for each i and g 6 SU(2). The gauge 
transformation G x . y associated with the translation also 
take the above form 



G x (i) = r\ x (i)g x 

Gyii) = vy(i)g y 



(C3) 



G x {i) = g x , 

G y {i) = {-T x 9y, 

Gt* (i) = V % t 9T*, 

Gpji) = (-y^g Pxy , 

G Pxg (i) = {~r H v l pxy gp, g 



(Cll) 



where g x ,9y,9T*,9P* y i9P X y e SU(2),n t = ±1, and 
i] pxy — ±1. Moreover, only rjt = — 1 and r\ vxy = — 1 lead 
to non vanishing ansatz. So there are only two kinds of 
SU{2) spin liquids, 
SU2A spin liquids: 



SU2B spin liquids 

^i,?'+m ( ) 



u^t , m = odd, 



(C12) 



(C13) 
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